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CROSS REFERENCES TO RELATED APPLICATIONS 
[0001] This application claims priority from United States Provisional Patent Application 
Ser. No. 60/229,134 filed on August 30, 2000. 



5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0002] The invention is related generally to the field of interpretation of measurements 
made by well logging instruments for the purpose of determining the fluid content and 
permeability of earth formations. More specifically, the invention is related to methods 
1 0 for using Nuclear, Resistivity and Nuclear Magnetic Resonance (NMR) measurements 
and /or measurements made with a formation testing tool or pressure tests made in 
laminated reservoirs for detenriining a distribution of sands, shales and fluids in the 
reservoir and estimating permeability of the reservoir. 

15 2. Background of the Art 

[0003] A significant number of hydrocarbon reservoirs include deep water turbidite 
deposits that consist of thin bedded, laminated sands and shales. A common method for 
evaluating the hydrocarbon content of reservoirs is the use of resistivity measurements. 
20 In interpretation techniques known in the art, typically one or more types of porosity- 
related measurement will be combined with measurements of the electrical resistivity (or 
its inverse, electrical conductivity) of the earth formations to infer the fluid content within 
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the pore spaces of the earth formations. The fractional volumes of connate water and 
hydrocarbons can be inferred from empirical relationships of formation resistivity Rt 
with respect to porosity and connate water resistivity such as, for example, the well 
known Archie relationship. In the Archie relationship fractional volume of water in the 
pore space is represented, as shown in the following expression, by Sw - known as "water 
saturation": 



where a and m are empirically determined factors which relate the porosity (represented 
10 by <D) to the resistivity of the porous rock formation when it is completely water-saturated 
(Ro), R w represents the resistivity of the connate water disposed in the pore spaces of the 
formation, and m represents an empirically determined "cementation" exponent, n is the 
saturation exponent. 

1 5 [0004] Relationships such as the Archie formula shown in equation (1) do not work very 

well when the particular earth formation being analyzed includes some amount of 
extremely fine-grained, clay mineral-based components known in the art as "shale". 
Shale typically occurs, among other ways, in earth formations as "dispersed" shale, where 
particles of clay minerals occupy some of the pore spaces in the hydrocarbon-bearing 

20 earth formations, or as laminations (layers) of clay mineral-based rock interleaved with 




(1) 
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layers of reservoir-type rock in a particular earth formation. 



[0005] In the case of dispersed shale, various empirically derived relationships have been 
developed to calculate the fractional volume of pore space which is capable of containing 
5 movable (producible) hydrocarbons. The fractional volume of such formations which is 
occupied by dispersed shale can be estimated using such well logging devices as natural 
gamma ray radiation detectors. See for example, M. H. Waxman et al, "Electrical 
Conductivities in Oil Bearing Shaly Sands", SPE Journal, vol. 8, no. 2, Society of 
Petroleum Engineers, Richardson, TX (1968). 

0 

[0006] In the case of laminated shale, the layers sometimes are thick enough to be within 
the vertical resolution of, and therefore are determinable by, well logging instruments 
such as a natural gamma ray detector. In these cases, the shale layers are determined not 
to be reservoir rock formation and are generally ignored for purposes of determining 
5 hydrocarbon content of the particular earth formation. A problem in laminated shale 
reservoirs is where the shale laminations are not thick enough to be folly determined 
using gamma ray detectors and are not thick enough to have their electrical resistivity 
accurately determined by electrical resistivity measuring devices known in the art. 

I [0007] Sands that have high hydrocarbon saturation are typically more resistive than 
shales. In reservoirs consisting of thin laminations of sands and shales, conventional 
induction logging tools greatly underestimate the resistivity of the reservoir: the currents 
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induced in the formation by the logging tool flow preferentially through the conductive 
shale laminations creating a bias towards a higher formation conductivity. This could 
lead to an underestimation of hydrocarbon reserves. 

5 [0008] One method for estimating hydrocarbon content of earth formations where shale 

laminations are present was developed by Poupon. See A. Poupon et al, "A Contribution 
to Electrical Log Interpretation in Shaly Sands", Transactions AIME, Vol. 201, pp. 138- 
145 (1959). Generally the Poupon relationship assumes that the shale layers affect the 
overall electrical conductivity of the earth formation being analyzed in proportion to the 

1 0 fractional volume of the shale layers within the particular earth formation being analyzed. 
The fractional volume is typically represented by V sh (shale "volume"). Poupon's model 
also assumes that the electrical conductivity measured by the well logging instrument will 
include proportional effects of the shale layers, leaving the remainder of the measured 
electrical conductivity as originating in the "clean" (non-shale bearing) reservoir rock 

1 5 layers as shown in the following expression: 

where R, represents the electrical resistivity (inverse of conductivity) in the reservoir rock 
layers of the formation and represents the resistivity in the shale layers. 

20 
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[00091 The analysis by Poupon overlooks the effect of anisotropy in the resistivity of a 
reservoir including thinly laminated sands and shales. Use of improper evaluation models 
in many cases may result in an underestimation of reservoir producibility and 
hydrocarbon reserves by 40% or more as noted by van den Berg and Sandor. Analysis of 
well logging instrument measurements for determining the fluid content of possible 
hydrocarbon reservoirs includes calculating the fractional volume of pore space 
("porosity") and calculating the fractional volumes within the pore spaces of both 
hydrocarbons and connate water. As noted above, Archie's relationship may be used. 

[0010] In thinly laminated reservoirs where the wavelength of the interrogating 
electromagnetic wave is greater than the thickness of the individual layers, the reservoir 
exhibits an anisotropy in the resistivity. This anisotropy may be detected by using a 
logging tool that has, in addition to the usual transmitter coil and receiver coil aligned 
along with the axis of the borehole, a receiver or a transmitter coil aligned at an angle to 
the borehole axis. Such devices have been well described in the past for dip 
determination. See, for example, United States patent 3,510,757 to Huston and United 
States patent 5,1 15,198 to Gianzero, 

[0011] Co-pending United States Patent Application Ser. No. 09/474,049 (the '049 
application) filed on December 28, 1999 and the contents of which are fully incorporated 
herein by reference, disclosed a method of accounting for the distribution of shale in a 
reservoir including laminated shaly sands using vertical and horizontal conductivities 



derived from multi-component induction data. Data such as from a borehole resistivity 
imaging tool give measurements of the dip angle of the reservoir, and the resistivity and 
thickness of the layers on a fine scale. The measurements made by the borehole 
resistivity imaging tool are calibrated with the data from the induction logging tool that 
gives measurements having a lower resolution than the borehole resistivity imaging tool. 
A tensor petrophysical model determines the laminar shale volume and laminar sand 
conductivity from vertical and horizontal conductivities derived from the log data. The 
volume of dispersed shale and the total and effective porosities of the laminar sand 
fraction are determined using a Thomas-Stieber-Juhasz approach. Removal of laminar 
shale conductivity and porosity effects reduces the laminated shaly sand problem to a 
single dispersed shaly sand model to which the Waxman-Smits equation can be applied. 

[0012] Co-pending United States Patent Application Ser. No. 09/539,053 (the '053 
application) filed on March 30, 2000, having the same assignee as the present application, 
and the contents of which are fully incorporated herein by reference, discloses a method 
of accounting for the distribution of shale and water in a reservoir including laminated 
shaly sands using vertical and horizontal conductivities derived from multi-component 
induction data. Along with an induction logging tool, data may also be acquired using a 
borehole resistivity imaging tool. The data from the borehole resistivity imaging tool 
give measurements of the dip angle of the reservoir, and the resistivity and thickness of 
the layers on a fine scale. The measurements made by the borehole resistivity imaging 
tool are calibrated with the data from the induction logging tool that gives measurements 



having a lower resolution than the borehole resistivity imaging tool. The measurements 
made by the borehole resistivity imaging tool can be used to give an estimate of V , , 
the volume fraction of laminar shale. A tensor petrophysical model determines the 
laminar shale volume and laminar sand conductivity from vertical and horizontal 
conductivities derived from the log data. The volume of dispersed shale, the total and 
effective porosities of the laminar sand fraction as well as the effects of clay-bound water 
in the formation are determined. 

[0013] The method of the '053 application is not readily applicable to reservoirs in which 
the sands may be intrinsically anisotropic without making additional assumptions about 
the sand properties. Sands in turbidite deposits commonly comprise thin laminae having 
different grains size and/or sorting: the individual laminae may be isotropic but on a 
macroscopic scale relevant to logging applications, the laminations exhibit transverse 
isotropy. In addition, a reservoir including turbiditic sands exhibits an anisotropic 
permeability. Being able to determine this anisotropic permeability is important from the 
standpoint of reservoir development. This is an issue not addressed in the '053 
application and of considerable importance in development of hydrocarbon reservoirs. 

SUMMARY OF THE INVENTION 
[0014] In one aspect of the invention, a method of petrophysical evaluation of a formation 
is disclosed wherein horizontal and vertical resistivities of the formation are inverted 
using a tensor petrophysical model to give a first estimate of fractional volume of 
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laminated shale in the formation. This first estimate of fractional volume of laminated 
shale is comparted to a second estimate obtained from measurements of density and/or 
neutron porosity of the formation using a volumetric model. If the second estimate of 
fractional shale volume is greater than the first estimate of fractional shale volume, the 
5 horizontal and vertical resistivities are inverted using a tensor petrophysical model 

including the second estimate of fractional shale volume and obtaining a vertical and 
horizontal resistivity of an anisotropic sand component of the formation. This vertical 
and horizontal resistivity of the anisotropic sand component is used in conjunction with at 
least one additional measurement selected from the group consisting of: of (i) NMR 
10 measurements of the formation, and, (ii) a bulk permeability of the sand component to 
obtain properties a coarse and a fine grain portion of the sand component. The obtained 
properties of the coarse and fine grain portions of the sand include water saturations, and 
resistivities. 

15 [0015] The properties of the coarse and fine grain portions of the sand are derived using 
an iterative solution process wherein a out of a family of possible distributions of said 
properties, a selection is made that matches the NMR measurement or the bulk 
permeability measurement. Relationships such as the Timur Coates equation may be 
used for the purpose. The bulk permeability measurement may be obtained from a 

20 formation testing instrument, a pressure build up test, a pressure drawdown test or from 
an NMR diffusion measurement. 
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[0016| Measurements of the horizontal and vertical resistivity may be obtained using a 
transverse induction logging tool, or from a conventional induction logging tool and a 
focused current resistivity tool. 

BRIEF DESCRIPTION OF THE FIGURES 
[0017] FIG. 1 (PRIOR ART) shows a resistivity imaging tool suspended in a borehole; 
FIG. 2 (PRIOR ART) is a mechanical schematic view of the imaging tool of Fig. 1; 
FIG. 2A (PRIOR ART) is a detail view of an electrode pad for the tool of Figs. 1, 2;; 
FIG. 3 (PRIOR ART) is a pictorial view of a composite imaging log obtained by 
merging the resistivity image data shown in acoustic image data; 
FIGS. 4A and 4B illustrate the principal steps of the process of one embodiment of the 
invention. 

FIG. 5Ais a schematic illustration of the distribution of water and hydrocarbons in a 
porous reservoir. 

FIG. SB shows the relationship between water saturation and the volume of coarse grain 
sand in a porous rock. 

FIG. 6 A is a plot of resisitivity index as a function of water saturation for a bimodal sand 
reservoir. 

FIG. 6B shows the effect of water saturation on the resistivity anisotropy of a bimodal 
sand. 

FIG. 7 A shows possible solutions for the inversion of measured resistivity of a bimodal 
sand reservoir. 

FIG. 7B shows the dependency of spherical permeability of a bimodal sand reservoir on 
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the volume fraction of the coarse grained component. 

DETAILED DESCRIPTION OF THE INVENTION 



[0018] The present invention is best understood by referring to Figs. 1 - 7. Fig. 4 is a 
schematic flowchart of the major steps of the process used in the present invention. 
Fig. 1 shows a composite imaging tool 10 suspended in a borehole 12, that penetrates 
earth formations such as 13, from a suitable cable 14 that passes over a sheave 16 
mounted on drilling rig 18. By industry standard, the cable 14 includes a stress member 
and seven conductors for transmitting commands to the tool and for receiving data back 
from the tool as well as power for the tool The tool 10 is raised and lowered by 
draw works 20. Electronic module 22, on the surface 23, transmits the required operating 
commands downhole and in return, receives digital data back which may be recorded on 
an archival storage medium of any desired type for concurrent or later processing. A data 
processor 24, such as a suitable computer, may be provided for performing data analysis 
in the field in real time or the recorded data may be sent to a processing center or both for 
post processing of the data. 

[0019] Fig. 2 is a schematic external view of the unified borehole sidewall imager 
system. This may be used to provide the data that may be used in an optional 
embodiment of the invention. The tool 10 comprising the imager system includes four 
important components: 1) resistivity arrays 26; 2) electronics modules 28 and 38; 3) a 
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mud cell 30; and 4) a circumferential acoustic televiewer 32. All of the components are 
mounted on a mandrel 34 in a conventional well-known manner. The outer diameter of 
the assembly is about 5.4 inches and about five feet long. An orientation module 36 
including a magnetometer and an inertial guidance system is mounted above the imaging 
assemblies 26 and 32. The upper portion 38 of the tool 10 contains a telemetry module 
for sampling, digitizing and transmission of the data samples from the various 
components uphole to surface electronics 22 in a conventional manner. Preferably the 
acoustic data are digitized although in an alternate arrangement, the data may be retained 
in analog form for transmission to the surface where it is later digitized by surface 
electronics 22. 

[0020] Also shown in Fig. 2 are three resistivity arrays 26 (a fourth array is hidden in this 
view). Referring to Figs. 2 and 2A, each array includes 32 electrodes or buttons identified 
as 39 that are mounted on a pad such as 40 in four rows of eight electrodes each. Because 
of design considerations, the respective rows preferably are staggered as shown, to 
improve the spatial resolution. For reasons of clarity, less than eight buttons are shown in 
Fig. 2A. For a 5.375" diameter assembly, each pad can be no more than about 4.0 inches 
wide. The pads are secured to extendable arms such as 42. Hydraulic or spring-loaded 
caliper-arm actuators (not shown) of any well-known type extend the pads and their 
electrodes against the borehole sidewall for resistivity measurements. In addition, the 
extendable caliper arms 42 provide the actual measurement of the borehole diameter as is 
well known in the art. Using time-division multiplexing, the voltage drop and current 
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How is measured between a common electrode on the tool and the respective electrodes 
on each array to furnish a measure of the resistivity of the sidewall (or its inverse, 
conductivity) as a function of azimuth. 

[0021] The acoustic imager that forms the circumferential borehole imaging system 32 
provides 360 ° sampling of the sidewall acoustic reflectivity data from which a 
continuous acoustic imaging log or sonogram can be constructed to provide a display of 
the imaged data. 

[0022] The borehole resistivity imaging tool arrays necessarily allow sampling only 
across preselected angular segments of the borehole sidewall. From those data, a 
resistivity imaging log, consisting of data strips, one strip per array, separated by gaps, 
can be constructed and displayed. The angular width of each data-scan strip is equal to 2 
sin {S/(2R)}, where S is the array width and R is the borehole radius. The common 
data from the two imagers are merged together in a data processing operation to provide a 
substantially seamless display as shown in Fig. 3. The merging incorporates equalizing 
the dynamic range of the resistivity measurements with respect to the acoustic 
measurements. That balance is essential in order that the continuity of a displayed textural 
feature is not distorted when scanning across a resistivity segment of the display, between 
adjacent acoustic segments. 

[0023] The display in Fig. 3 incorporates measurements from directional sensors to align 
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the resistivity measurements with geographical coordinates (North, East, South, West), 
with the resistivity image being "unfolded" to provide a flat image of the cylindrical 
surface of the borehole. Those versed in the art would recognize that when a plane 
intersects a circular cylinder at an angle, the unrolled image of the plane would appears as 
a sinusoid. The display in Fig. 3 shows many such sinusoids, some corresponding to 
bedding planes and others corresponding to fractures. The dip angle and the dip direction 
corresponding to the various sinusoids are determined in the present invention using 
known methods. When these data are combined with measurements from other logs, 
such as a gamma ray or a neutron log, discrete layers of different lithologies may be 
identified. In particular, over a gross interval of the order of several meters or so, the 
fractional volume of laminated shale present in a laminated reservoir may be determined. 

[0024] With flat dips, the sinusoids have essentially zero amplitude. In one aspect of the 
present invention, the resistivity measurements are averaged circumferentially and 
vertically within each identified layer to give an average resistivity measurement for each 
layer identified above. Once this is done, the subsurface may be characterized by a 
number of plane layers, each of which has a constant resistivity. With the resolution of 
the button-electrode tool, these layers may range in thickness from a few millimeters to a 
few centimeters. 

[0025] Those versed in the art would recognize that when the bed boundaries are dipping, 
then the currents into the electrodes, particularly those in the dip direction, on the pads 
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may not be confined to a single layer and hence not represent the resistivity of the layer at 
the borehole. In one aspect of the invention, the averaging described above is limited to 
electrodes in the strike direction: these measurements would be more likely representative 
of the true formation resistivity at the depth of measurement. 

[0026] The resistivity measurements obtained by the averaging process correspond to 
layers that are beyond the resolution of electromagnetic induction logging tools or 
propagation resistivity tools. Accordingly, the resistivity measurements obtained at this 
point are averaged to give resistivities on a scale that would be measurable by an 
induction logging tool. 

[0027] As would be known to those versed in the art, a finely laminated sequence of 
layers having different resistivities exhibits a transverse isotropy on a larger scale where 
the wavelength of the electromagnetic wave is much greater than the layer thickness. 
This condition is easily satisfied even for propagation resistivity tools that, e.g., operate 
at a frequency of 2MHz (with a wavelength X » 6 meters); for induction logging tools 
that have frequencies of the order of 50kHz to 200kHz, the wavelengths are even longer. 
For such interrogating frequencies, the layered medium is characterized by a horizontal 
resistivity R h m and a vertical resistivity R v * given by: 




(3) 
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and 



KJ (4) 

where ^ is a window used to average the resistivities, Ah is the depth sampling interval 
of the electrodes, and R t is the measured resistivity for a given depth. 

(0028] In this invention, the terms "horizontal" and "vertical" are to be understood in 
terms of reference to the bedding planes and the anisotropy axes of the subsurface 
formations, i.e., "horizontal" refers to parallel to the bedding plane, and "vertical" refers 
to vertical to the bedding plane. Where the beds of the formation are dipping, the 
anisotropy axis is taken to be the normal to the bedding plane. When the borehole is 
inclined to the bedding plane, data from the orientation module 36 in Fig. 1, may be 
used to correct the resistivity measurements made by the resistivity imaging tool to give 
measurements parallel to and perpendicular to the bedding planes. 

[0029] Those versed in the art would recognize that the resistivity measurements made by 
the electrode-pad system described above may be in error and, in particular, may need to 
have a scaling factor applied to the data. When this data is acquired, it may be calibrated 
by relating the values given by equations (3) and (4) to data from an induction logging 
tool or a propagation resistivity tool. 
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(0030| Referring now to Fig. 4, one optional embodiment of the invention starts with 
data acquired by a borehole resistivity imaging tool such as is described in United States 
Patent 5,502,686 issued to Dory et al., and the contents of which are fully incorporated 
here by reference. It should be noted that the Dory patent is an example of a device that 
can be used for obtaining measurements borehole resistivity measurements: any other 
suitable device could also be used. The process of the invention starts with an initial 
model 101 for the structure of the reservoir. This initial model comprises a laminated 
shale fraction and a sand fraction. This initial model may be derived from the resistivity 
imaging tool described above. A horizontal and vertical conductivity C sh . lamh and 
c sh-iam.v ) of the shale fraction is assumed or is measured 103; if measurements are to be 
made within a borehole, this may be done by using a Transverse Induction Logging Tool 
(TILT) on a thick section of shale in proximity to the reservoir. The resistivity of the 
"bulk" shale may also be obtained from core measurements. 

[0031] An induction or wave propagation tool is used to make measurements of the 
vertical and horizontal resistivity of the earth formations. For example, United States 
Patent 5,781,436 to Forgangetal , the contents of which are fully incorporated here by 
reference, discloses a method an apparatus for making measurements of horizontal and 
vertical resistivities of a transversely isotropic formation. 

[0032] The method disclosed by Forgang et al comprises selectively passing an 
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alternating current through transmitter coils inserted into the welibore. Each of the 
transmitter coils has a magnetic moment direction different from the magnetic moment 
direction of the other ones of the transmitter coils. The alternating current includes a first 
and a second frequency. The amplitude at the first frequency has a predetermined 
relationship to the amplitude at the second frequency. The relationship corresponds to the 
first and the second frequencies. The method includes selectively receiving voltages 
induced in a receiver coil having a sensitive direction substantially parallel to the axis of 
the corresponding transmitter coil through which the alternating current is passed. A 
difference in magnitudes between a component of the received voltage at the first 
frequency and a component of the voltage at the second frequency is measured, and 
conductivity is calculated from the difference in magnitudes of the components of the 
received voltage at the two frequencies. The Forgang patent is cited only by way of 
example of an induction device for obtaining horizontal and vertical resistivities of a 
formation and there are other teachings on obtaining these properties of subterranean 
formation. 

[0033] An example of a propagation resistivity tool for making measurements of 
horizontal and vertical resistivities is described by Rosihal (U.S. Patent 5,329,448) 
discloses a method for determining the horizontal and vertical conductivities from a 
propagation logging device. The method assumes that 0, the angle between the borehole 
axis and the normal to the bedding plane, is known. Conductivity estimates are obtained 
by two methods. The first method measures the attenuation of the amplitude of the 
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received signal between two receivers and derives a first estimate of conductivity from 
this attenuation. The second method measures the phase difference between the received 



shift. Two estimates are used to give the starting estimate of a conductivity model and 
based on this model, an attenuation and a phase shift for the two receivers are calculated. 
An iterative scheme is then used to update the initial conductivity model until a good 
match is obtained between the model output and the actual measured attenuation and 
phase shift. 

[0034] As described in the '053 and '967 applications, measurements R lh and R l v made by 
TILT or other suitable device are inverted 111 to give an estimate of the laminar shale 
volume and the sand conductivity, assuming that the sand component is isotropic. In 
terms of resistivity, 



signals at two receivers and derives a second estimate of conductivity from this phase 



K + *„.,..)+ (rS - **.,.,) • V(T^)} (5) 



where 




AR = 4-R, 



tiso 
sd 




(6) 



R Jrf &0 is the 'isotropic' sand resistivity. If the shale is isotropic, ( R shh = R skv ), then this 
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resistivity is identical to the sand resisitivity. AR is the correction for anisotropic shale. 
AR becomes zero for an isotropic shale ( R sh h = R xh v ). 



[0035] The inversion also gives 



(7) 



An independent estimate of the laminar shale volume V sM ^ is obtained 107 from 
volumetric measurements using density or neutron logs 105 and using a method such as 
the well-known Waxman-Smits or Thomas-Stieber methods. Acoustic imaging logs may 
also be used to get the volume fraction of laminar shale. Obtaining this independent 
estimate of laminar shale volume V sMt ^ would be known to those versed in the art and is 
not discussed further here. 



[0036] Next, a check is made 113 to see if the V sU from the TILT inversion agrees with 
the V sh4t from neutron or density logs. If the two estimates of laminar shale volume are 
consistent 115, then the assumption of an isotropic sand at 111 is valid and classical 
methods based on Archie or Waxman-Smits are used to determine water saturation of the 
sands. Alternatively, the Dual- Water method disclosed in the '053 application are used to 
determine the water saturation in the sand. 
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[0037| If the answer at 1 13 is "No", then a check is made to see if V ihA ra is greater than 
V stu 117. If the answer is "Yes", then this is an indication to change 1 19 the 
assumptions made at 103. After changing the assumptions about the horizontal and 
vertical shale conductivities, the process goes back to 101. If the vertical shale 
conductivity is obtained from an actual measurement, then the assumed shale anisotropy 
factor is in error and a new value is chosen. Alternatively, the input parameters for the 
Thomas-Stieber calculations must be modified. 

[0038] If the answer at 117 is "No", then this is an indication that the sands component is 
anisotropic 121. In this case, the TILT resistivity data are inverted 123 using the value of 
Vsm. ts obtained at 107, e.g., using Thomas-Stieber and the method of the '049 application 
or the method of the '053 application, to give a water saturation S w . 

[0039] As a check on the determination of water saturation from the inversion of TILT 
data, NMR data are obtained 127 and from the NMR data, a determination of the bulk 
volume of irreducible water in the formation is made. Methods for determination of 
irreducible water saturation S wjrr from NMR data are disclosed in US patents 5,412,320 
and 5,557,200, the contents of which are fully incorporated herein by reference and not 
discussed further here. The bulk irreducible water content is given by S wirr (f> H, where <p 
is the porosity and H is the thickness, while the bulk water content from the TILT data is 
given by S w <pH. A check is made 129 to see if the former quantity is less than the latter 
quantity. If the test at 129 is negative, then it is an indication that there is a problem in 
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the TILT model and the process is restarted at 103. If the test at 129 is positive, then use 
is made of horizontal and vertical sand resistivities R sdh and R sd v in the TILT inversion 
123 given by 



(8) 



R sd,v = 



R 



/,v 



1 Vsh-lJS 



(9) 



[0040] The sand anisotropy with resistivity values R sdh and R sd v is indicative of a 
laminated sand layer. These values of R sdJt and R sd v are inverted to give a layered model 
131 comprising isotropic sand layers and the laminated shale component determined 
above. In order to perform this inversion, an estimate of the number and thicknesses of 
the sand layers is required. This may be obtained from a resistivity imaging tool as 
discussed in the '053 and '967 applications or it may be obtained using NMR data 127. 
From the distribution of relaxation times T, and T 2 of NMR data, a distribution of volume 
fractions of individual sand components 133 may be obtained using known methods. 
Alternatively, core information or sedimentologic information about the reservoir may be 
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used to give the volume fractions of the sand components. 

[0041 1 For the particular case of two components, eqns. (8) and (9) may be inverted to 
give a unique solution Tor isotropic conductivities for two sand layers having the 
5 appropriate volume fractions. If the number of sand layers in the laminated sand 

component is more than two, then there is no unique model of isotropic sand layers 
having an observed vertical and horizontal resistivity on a macroscopic scale. For the 
case of more than two sand layers, additional information, such as some relationship 
between the individual sand resistivities, is necessary to obtain the laminated sand 
10 component of the reservoir. The result of this inversion 131 of a model of volume 

fractions (V sh . lam , V*,, . . ) and resistivities ( , v , m h , R, ^ , , R u sd 2 . .). 
For convenience, the discussion hereafter is limited to two sand layers though it is to be 
understood that additional layering of the laminated sand component is possible. 

1 5 [0042] Using assumed values for the water saturated sand in horizontal and vertical 

direction Ro^ h and Ro tSiv the water saturation of the individual sand layer S WJ (i-th layer) 
are calculated separately 135 using the layer resistivity R^ obtained at 131. Depending 
on the saturation equation (Archie-equation, Waxman-Smits- equation) the following 
parameters are necessaiy as input: (i) Formation water resistivity, (ii) porosity or 

20 formation factor of the layer, (iii) saturation exponent of the layer, and, (iv) Waxman- 
Smits-parameters in case of dispersed shale in the sand layer. 
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[0043| A direct method for saturation calculation for a bimodal sand is given in Schoen et 
al. As shown therein, a simple model is presented to describe the influence of changing 
water content upon the electrical properties for a thinly laminated bimodal Archie-type 
sand. The sand comprises a coarse-grained component and a fine-grained component. The 
5 coarse sand fraction is characterized by low irreducible water saturation and the fine 
fraction, by high irreducible water saturation. As shown in Fig, 5a, a single sand layer 
may be depicted by four constituents: the matrix 201, capillary bound water 203, movable 
water 205 and hydrocarbon 207. The volume fraction of 203, 205 and 207 is the porosity 
of the sand. The change of water content only takes place in the pore space occupied by 

1 0 the mobile fluid pore space (the combination of 205 and 207). In this example, the pore 
space is occupied by immobile water (capillary bound water) and movable fluids 
(movable water, oil, and gas). The complete pore space therefore is not available for the 
variation of the water saturation, but only the mobile water fraction, 1- S w>irT . Detailed 
information about the saturation behavior in the individual layers could be derived from 

15 capillary pressure curves. For the simplified forward calculation this variation is 

expressed by a parameter defined as p and the total water saturation, assuming no clay- 
bound water, component can be written: 

S w = S wjrr + p • (l - S WM ) = /? + S w>irr (1-0) (10) 

20 

The parameter P describes how the "free water pore space" is occupied; for P=0, S w = S 



24 



while for P = 1, S w = /. 



[0044] For determining a variation of the total water content, one embodiment of the 
present invention assumes that all sand fractions are characterized by the same value of 
for P but different S Wttrr Thus, the "free water pore space" is occupied in the same ratio 
for all sand fractions as it is described schematically in Fig. 5b. Shown in Fig. 5b is a 
plot of the water saturation (ordinate) as a function of the percentage of coarse grained 
component of the sand (abscissa). Five different values of P are shown in the curves 
221a, 221b. . . 22 le, the last of which corresponds to p = 0 and having a value of 

S W = $ w.irr 



[0045] For a two-component sand interval (e.g., coarse, subscript c and fine, subscript f) 
we have two different values for S w irr (S wjrrx and S WJrrJ ) and the vertical and horizontal 
resistivities are 



(11) 



and 
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R t.h - R w 



(12) 



where V a fry- are the volume fractions of the coarse and fine sand, F c v , F c h F f v , F f h are the 
formation factors for coarse and fine sand related to vertical and horizontal current 
direction. 

5 

[0046] The mean saturation of the interval is 



1 0 where <p C9 <p f are the porosities of the coarse and fine sand fraction. In a preferred 

embodiment of the invention, it is assumed that the two sand layers have the same 
intrinsic microscopic anisotropy of the formation factor. In most cases, this microscopic 
anisotropy of the individual sand layers can be neglected when compared with the 
macroscopic anisotropy effect. Then the resistivity horizontal and vertical resistivity 
15 indices are 

26 



Rh = iK • IP + s wirrc (i- /?)] + V/[P+ s wjrrJ • (i- /?)]"}"' ( i4) 

and 

v \fi + ^ • (1 - ;W + + S WMJ • (1 - /?)]" (15) 



[0047] The different saturations in the two sand fractions is electrically described now by 
the two resistivity indices. With Eq. (13) the resistivity versus water saturation 
relationships can be calculated. Fig. 6a shows a logarithmic plot of the calculated 
resistivity index (ordinate) versus water saturation (abscissa). For each individual sand 
layer, a saturation exponent index of n = equal to 2 was assumed. For the Archie 
relationship, the curve is a straight line 241 with a slope of 2. 



[0048] All the curves in Fig. 6a correspond to an S Wt irr c =0. 1 . The curves 243a, 243b are 
the vertical and horizontal resistivity indices for a value of S Wt irr/ = 0.2; the curves 245a, 
245b are the vertical and horizontal resistivity indices for a value of S w , rr/ = 0.3; while the 
curves 247a, 247b are the vertical and horizontal resistivity indices for a value of S w trrJ = 
0.5. 



[0049] The anisotropy effect associated with the different partial water saturations 
however results in the following: 

(1) The saturation exponent relationship is a curve not a straight line and the 
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saturation index exponent, «, of for the "composite" is dependent on water 
saturation. 

(2) The vertical resistivity indices are much higher, and the horizontal resistivity 
index indices are smaller when compared to values for an Archie sand with n = 2. 

(3) The deviation of the resistivity index saturation exponent from n = 2 increases 
with decreasing water saturation (increasing hydrocarbon saturation). 

(4) For S w equal to 1, both curves converge to the Archie n = 2 curve. 

(5) All the saturation exponent variations are controlled by the contrast in irreducible 
water saturation and volume fraction of the sands. 

[0050] The ratios of resistivity anisotropy versus water saturation for the case shown in 
Fig. 6a are shown in Fig. 6b for the same cases shown in Fig. 6a. As in Fig. 6a, S w irr c 
=0. 1 . The curves 261, 263 and 265 correspond to values of S w> irrJ equal to 0.2, 0.3 and 
0.5 respectively. The abscissa is the mean water saturation while the ordinate is the 
resistivity anisotropy ratio. Fig. 6b illustrates the remarkable influence of low water 
saturation (high hydrocarbon saturation) on anisotropy. The maximum possible 
anisotropy ratio is determined by the maximum possible contrast of the two irreducible 
water saturations. Jn this example, maximum anisotropy strongly trends to the higher 
"mean" water saturation as a result of the high value S w irrf 

[0051] Having reviewed the effect of water saturation on the anisotropy of a laminated 
sand, we now address the problem of inverting measured values of horizontal and vertical 
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resistivities for a laminated sand to obtain water saturations. The saturation 
determination problem for laminated, bimodal sands can be described as follows: 

( 1 ) For the laminated macroscopic anisotropic sand, no constant saturation exponent 
n can be defined. 

(2) A formalistic application of Archie's equation with a constant saturation exponent 
results in different saturation values for horizontal and vertical measured 
resistivities, and a physically incorrect result. 

[0052] In a preferred embodiment of the invention, a hydrocarbon-bearing bimodal sand, 
is modeled as comprising an alternating square step profile of a coarse (subscript c) and a 
fine fraction (subscript f). Other embodiments of the invention use other profiles, such as 
a triangular gradation of grain size between two limits. The calculations for such profiles 
would be readily obtainable for those versed in the art based upon discussions given here 
and in Schoen, Mollison & Georgi (1999) and are not discussed further. 

[0053] The pore space of both sand layers has different water saturation. The water 
saturation of the fine sand S w/ is higher than the water saturation of the coarse sand S w c . 
In particular^different grain sizes result jn proportionately different pore sizes and can be 
directly related to irreducible water saturation and capillary pressure. For both sands we 
assume that Archie's law is valid and no shale is present. The mean water saturation of 
the interval S w is: 

S^K'S^ c + V r S wJ (16) 
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The horizontal and vertical resistivities of the composite sediment are determined as 
follows 



1-1 



' S wJ 



(17) 



S" 



08) 



where F c vt F c h F /Vf F f h are the formation factors for the coarse and fine sands for the 
vertical and horizontal current directions. 

[0054] The following equations focus on the influence of the water saturation upon the 
10 resistivity and anisotropy. For well-sorted sands the porosity is independent of grain size; 
thus, a reasonable assumption is that the fine and coarse sand fraction porosities are 
similar. Further, it is attractive to assume that the formation factor is a simple single 
valued scalar. Therefore, in a first approximation, in the preferred embodiment of the 
invention, we assume that within the interval the formation factor is constant within the 
1 5 interval of interest. This assumption results in a simplification of the equations. 

[0055] For the bimodal sand we derive three equations: horizontal resistivity index (Eq. 
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19), vertical resistivity index (Eq. 20), and volumetric closure (Eq. 21): 




\v S" + 




(19) 



R, 




(20) 



(21) 



The inversion process results in two water saturation values, S w c and S w/ . This calculation 
requires the volume fraction of one component (e.g., coarse component V c or the fine 



[0056] Because n is not constant and Archie's law is not valid for the composite 
resistivity index, the following calculation is done to find one consistent solution for the 
water saturation values of the two layers. We start with eqs. (19) and (20). From our a 
priori information of the coarse component, S w c n from eq. (19) is inserted into eq. (20) , 
resulting in the water saturation of the fine sand fraction: 




RI1 [ + RI~ l (\-2V V 


+ 


I i-K J 


ll"*/; 1 + */;'• (i-2r,) 


- 2 


4 L !-^c 



(22) 



Note that only the positive root is physically realistic. 

The water saturation of the coarse fraction can then be determined as 

l 



w,c 



V rtn 



Ri,- S . 



»J J 



(23). 



Thus, both saturation values are obtained and the mean saturation is obtained from 
(16). 



[0057J Returning now to Fig. 4B, in an optional embodiment of the invention, the results 
derived above are used for forward modeling of the properties of the reservoir rock 139 
and a consistency check may be made going back to 103. Still referring to Fig. 4B, an 
anisotropic permeability calculation is made 141 for the reservoir with the properties as 
derived above. This is accomplished as described immediately below. 
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[0058| The starting point for the permeability determination is the Coates-Timur equation 
'*Y (<j>-BVl\ b 



k = 



BVI 



(24). 



An alternate form of the equation is given by Kenyon as Je= C<f> a T b and may be used. 

For the bimodal sand distribution derived above, for the two individual layers, the 
microscopic permeabilities are 



BVI 



f J 



(25) 



and 



, JiX U C -BVI\ 



(26) 



The macroscopic permeabilities in the horizontal and vertical directions are 

**■*;•*. + */•*/ (27) 



and 



(28). 
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Solution of cqs. (19) - (28) requires knowledge of the volume fractions of the coarse- and 
fine-grained components of the sand. This may be obtained from core analysis or from 
distribution of NMR relaxation times. 



[0059] In an alternate embodiment of the invention, eqs. (19)- (28) are solved for a suite 
of values of V c and V f to give a family of solutions. Each of these has associated values 
of k h and From these derived horizontal and vertical permeabilities, a global 
permeability is determined that averages over the individual layers and directions. One 
such global measure of permeability is the spherical permeability given by 



[0060] In a preferred embodiment of the invention, it is assumed that the permeability 
derived from NMR measurements is this spherical permeability. Alternatively, a 
permeability may be obtained using the formation testing instrument marketed by Baker 
Hughes under the mark RCI™ . Bulk permeability may also be obtained from reservoir 
pressure buildup tests of pressure drawdown tests. 

[0061] Having a measurement of the bulk permeability gives a unique solution to the 
inversion problem. In one embodiment of the invention, this solution is obtained using a 
nonlinear iterative algorithm. Such an algorithm would be known to those versed in the 
art and is not discussed further. Alternatively, the unique solution may be obtained by a 
table look up in a processor or by an equivalent graphical solution given in the example 




(29) 
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below. 



[0062] In order to test the method a data set was created by forward calculation. With the 
resulting data then the iterative inversion algorithm was started. 

Data for forward calculation: 
Rw = 0.1 Ohm m m = n = 2 

V c = 0.4 V f = 0.6 <D C = 0.34 O f =0.25 BVI C = 0.05 BVI, = 0.10 
using Coates equation the microscopic permeabilities are : k f = 88 md k c = 4495 md 

Result of forward calculation (macroscopic model data): 

R„ = 22.0 Ohm m R h = 14.3 Ohmm 

0 = 0.286 BVI = 0.08 

k h = 1851 md k v =145md k sph = 791md 

Result of forward calculation (macroscopic model data): 

R v = 22.0Ohmm R^ 14.3 Ohmm 

<D = 0.286 BVI = 0.08 

k h = 1851md k v =145md ^ = 791 md 

[0063] Fig. 7a shows the relationship between BVI C and BVI f for a range of assumptions 
of Kc between 0.1 and 0.6, i.e., all the solutions fit the measured resistivity values R, and 
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(0064] Turning now to Fig. 7b, three curves 311, 313 and 315 are shown of the spherical 
permeability (ordinate) as a function of Vc. The three curves all have the same average 
porosity of 0.286 but the porosity in the coarse fraction is 0.33 for 315, 0.34 for 313 and 
0.35 for 311. The range of values from 0.33 to 0.35 is a typical "noise" in the 
determination of porosity. A measured spherical porosity value is indicated by 321. It 
may be seen in Fig. 7b that a range of values between 325 and 327 for Vc are possible 
solutions to the inversion problem. Table I gives the corresponding reservoir properties. 
TABLE I: Reservo ir properties derived from combined analysis for three different 



porosity assumptions 





v e 


BVI C 


Mmd) 


kh(md) 


K 


*r 


V f 


BVI f 


Mmd) 


Mmd) 




0.33 


0.47 


0.052 


3318 


1614 


8.5 


0.26 


0.53 


0.104 


104 


191 




0.34 


0.40 


0.050 


4493 


1850 


13 


0.25 


0.60 


0.100 


88 


145 




0.35 


0.37 


0.046 


6619 


2039 


19 


0.24 


0.63 


0.095 


76 


109 ! 





(0065) It is worth noting that the determined permeability anisotropy (defined as a ratio 
of the horizontal to vertical permeability) ranges from 8.5 to 19. In comparison, the 
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resistivity anisotropy, defined as the ratio of the vertical to horizontal resistivity, is 
(22.0/14.3) - 1.54. Those versed in the art would recognize that a commonly used form 
of application of the Coates eqn. (26) is 



k = 



,™ J2 4-bviY 

10M 'bvTJ (30) 



For BVI « <|>, the result is 



16 

At = io 4 -— on 

BVI 2 { ; * 



When combined with empirical relationships of Klein between permeability and 
resistivity, the result is 

K = \0*.^-R h k h = lO*.^-.R v (32) 



This would imply that the resistivity and permeability anisotropy are equal, in 
contradiction to the results derived here. This apparent paradox results from the 
combination of macroscopic model equations for laminated materials with microscopic 
empirical correlations between hydraulic connectivity and electrical resistivity. 
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[0066] While specific embodiments of the microresistivity tool and induction logging 
tool have been discussed above, it is to be understood that the tools may be used either on 
a wireline or in an MWD environment. It is to be further understood that the anisotropy 
measurements discussed above with reference to an induction logging tool may also be 
obtained using a propagation resistivity tool. Specifically, in a preferred embodiment of 
the invention, a transverse induction logging tool has been described for obtaining 
measurements indicative of horizontal and vertical resistivities of formation. In an 
alternate embodiment of the invention, the horizontal resistivities may be obtained by a 
conventional induction logging tool with a coil axis parallel to the borehole axis while 
vertical resistivities may be obtained from measurements using a focused current logging 
tool. Such focused current logging tools would be known to those versed in the art and 
are not discussed further. 

[0067] While the foregoing disclosure is directed to the preferred embodiments of the 
invention, various modifications will be apparent to those skilled in the art. It is intended 
that all variations within the scope and spirit of the appended claims be embraced by the 
foregoing disclosure. 



38 



